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Abstract An optimized digital shaping filter has been 
developed for the Gerda experiment which searches 
for neutrinoless double beta decay in 76 Ge. The Gerda 
Phase I energy calibration data have been reprocessed 
and an average improvement of 0.3 keV in energy reso¬ 
lution (FWHM) at the 76 Ge Q value for Ozyd/3 decay is 
obtained. This is possible thanks to the enhanced low- 
frequency noise rejection of this Zero Area Cusp (ZAC) 
signal shaping filter. 

Keywords germanium detectors • enriched 76 Ge • 
neutrinoless double beta decay • signal processing • 

PACS 23.40.-s ft decay; double /3 decay; electron and 
muon capture • 14.60.St non-standard-model neutrinos, 
right-handed neutrinos, etc. • 29.40.Wk solid-state 
detectors • 29.85.-c computer data analysis 


1 Introduction 

Gerda (GERmanium Detector Array) [1] searches for 
neutrinoless double beta decay (Ozz/3/3 decay) in 76 Ge. 
The experiment is located at the underground Gran 
Sasso National Laboratory (LNGS) of INFN, Italy. Crys¬ 
tals made from isotopically modified germanium with 
a fraction of ~86% of ,6 Ge for a total mass of ~20 kg 
are operated as source and detector of the process. 

Several extensions of the Standard Model of particle 
physics predict the existence of Oi//3/3 decay, a process 
which violates lepton number conservation by two units 
and which is possible if neutrinos have a Majorana mass 
component. Ojz/ 3/3 decay is therefore of primary interest 
in the field of neutrino physics. Neglecting the nuclear 
recoil energy the energy released by a Ozz/3/3 event is 
shared by the two emitted electrons. Both electrons are 
stopped within ~1 mm of germanium and thus all avail¬ 
able energy is deposited in a small region inside the 
detector. Since distortions by bremsstrahlung are ex¬ 
pected to be small the 0 vpfi decay signature is a peak 
in the energy spectrum at the Q value of the reaction, 
Qpp, amounting to 2039 keV for 76 Ge. The most recent 
result of this process for 76 Ge was published by the 
Gerda collaboration with a 90 % confidence level (CL) 
limit on the Oiz/3/3 half-life of T®J 2 > 2.1 • 10 25 yr [2]. 

The sensitivity for detection of a possible Ovfifi de¬ 
cay signal depends on the total efficiency e (~75 % for 
Gerda Phase I), the enrichment fraction fe and the 
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isotopic mass mi of the considered isotope, the total 
source mass M, the background level and the energy 
resolution. The expected number of signal events n$ 
for a given half-life Ty 2 is [3]: 


n s 
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( 1 ) 


where N A is the Avogadro number and t the live time of 
the measurement. The expected number of background 
events ub within an energy window AE is: 


n B = BI ■ AE ■ M -t (2) 

with BI being the background index in cts/(keV-kg-yr). 
The size of AE is proportional to the energy resolu¬ 
tion at Qf 3 / 3 , expressed as full width at half-maximum 
(FWHM). The energy resolution is of primary impor¬ 
tance for the enhancement of the sensitivity and the 
modeling of background sources. If the event waveforms 
are fully digitized with enough band width, the opti¬ 
mization of energy resolution through a digital signal 
processing is possible. 

A new energy reconstruction shaping filter leading 
to an improved energy resolution has been developed 
(section 3), that is denoted as Zero Area Cusp (ZAC) 
filter. The Gerda experiment (section 2), the readout 
of the data (section 2.1) and the signal processing (sec¬ 
tion 2.2) are described first. After the optimization of 
the ZAC filter (section 4) the Phase I data have been 
reprocessed (section 5). 


2 The GERDA experiment 

The design and the construction of Gerda were tai¬ 
lored to background minimization. The germanium de¬ 
tectors are mounted in low-mass ultra-pure copper hold¬ 
ers and are directly inserted in 64 m 3 of liquid argon 
(LAr) acting as cooling medium and shield against ex¬ 
ternal background radiation. The argon cryostat is com¬ 
plemented by a water tank with 5 m diameter which 
further shields from neutron and gamma backgrounds. 
It is instrumented with photomultipliers to veto the cos¬ 
mic muons by detecting Cerenkov radiation. A further 
muon veto is provided by plastic scintillators installed 
on the top of the structure. A detailed description of 
the experimental setup is provided in Ref. [1], 

A first physics data collection, denoted as Phase I, 
was carried out between November 2011 and June 2013. 
In Phase I eight p-type senri-coaxial detectors enriched 
in 76 Ge from the Heidelberg-Moscow (HdM) [4] and 
IGEX [5] experiments and five Broad Energy Germa¬ 
nium (BEGe) detectors were used [6]. Three coaxial 
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Fig. 1 Typical readout scheme of a germanium detector 
with a charge sensitive preamplifier with open loop gain A. 
The detector with capacitance Cn is operated with inverse 
bias voltage HV. The charge Qi n is collected on the capaci¬ 
tor Cf which then discharges because of the presence of the 
feedback resistor Rf. 


detectors with natural isotopic abundance from the Ge¬ 
nius Test Facility (GTF) project [7,8] were also in¬ 
stalled. In a second physics run (Phase II) 30 BEGe 
detectors will be operated in addition to the eight semi¬ 
coaxial together with instrumentation to detect the LAr 
scintillation light to actively suppress background [9,10, 
11 ]. 

2.1 Signal readout and shaping with germanium 
detectors 

The typical readout of a germanium detector operated 
as a diode with inverse bias voltage applied consists 
of a charge sensitive preamplifier whose output wave 
form is either shaped and then processed by an analog 
to digital converter or, as in Gerda, directly digitized 
by a flash analog to digital converter (FADC). Fig. 1 
presents the detector and the charge sensitive pream¬ 
plifier system consisting of a junction gate field-effect 
transistor (JFET) coupled to a feedback circuit. The 
capacitor Cf integrates the charge from the detector 
causing a steep change in voltage at the preamplifier 
output. In order not to saturate the dynamic range of 
the preamplifier a feedback resistor Rf is connected in 
parallel to the capacitor to bring back the voltage to 
its baseline value. The shape of the preamplifier output 
pulse will then be characterized by a fast step, with rise 
time of about 0.5-1.5 ps corresponding to the charge 
collection process followed by an exponential decay with 
time constant r = RfCf. The values of Rf and Cf for 
the Gerda preamplifiers are 500 Ml? and ~0.3 pF, re¬ 
spectively, for a t of about 150 ps. A description of the 
Gerda readout scheme is given in Ref. [1]. 

Fig. 2 shows the signal and main intrinsic noise 
sources in the detector and preamplifier system. The in¬ 


noiseless 

preamplifier 



Fig. 2 Signal and main noise sources in a germanium detec¬ 
tor readout system. The trace recorded by the digitizer can be 
modeled as the output of a noiseless preamplifier, connected 
to a noiseless detector with capacitance Cn, a series volt¬ 
age generator and a parallel current generator with spectral 
densities s(/) and p(f), respectively. Q ■ S(t) is the original 
current signal and Ci is the preamplifier input capacitance. 


trinsic equivalent noise charge (ENC) for a given shap¬ 
ing time t s is given as: 

ENC 2 = a—C$+pA f C%+ 7 ( e(I G +I L ) + ^)T s 

9mT s \ Rf J 

( 3 ) 

where g m the JFET transconductance, k is the Boltz¬ 
mann constant, and T the operational temperature. 
The constants a, /3 and 7 are of order 1 depending on 
the signal shaping filter (c.f. Ref. [12]). The series noise 
(first term) is proportional to the total capacitance Ct 
which is the sum of the detector capacitance Cd , the 
feedback capacitance C / and the preamplifier input ca¬ 
pacitance Ci. The second term represents the 1// noise 
of the JFET with amplitude A f and is also proportional 
to the total capacitance. The third term is the parallel 
noise generated by the detector leakage current I L , the 
gate current Iq and the thermal noise of the feedback 
resistor Rf. The parallel noise is proportional to t s and 
the series noise to its inverse while the 1 // noise is in¬ 
dependent of t s . Therefore, the optimal shaping time 
is the one which minimizes the sum of the series and 
parallel noise. More detailed descriptions of the noise 
origin and its treatment in germanium detectors can 
be found in Refs. [12] and [13]. 

In Gerda Phase I an additional low-frequency dis¬ 
turbance comes from microphonics related to mechani¬ 
cal vibrations of the long contacts (30-60 cm) connect¬ 
ing the detectors to the preamplifiers. 

2.1.1 Digital Shaping 

In Gerda Phase I the signals were digitized with 14 bits 
precision and 100 MHz sampling frequency [1]. 16384 
samples were recorded per pulse (Fig. 3). After a ~80 ps 
long baseline the charge signal rises up with a ~1 ps 
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Fig. 3 Typical wave form recorded in Gerda Phase I. A 
~80 ps long baseline is recorded before each signal. The ex¬ 
ponential decay tail is from the discharge of the feedback 
capacitor. 


rise time followed by a ~80 |xs long exponential tail due 
to the discharge of the feedback capacitor. 

The energy estimation was performed by applying 
a shaping filter to the digitized signal. The advantages 
with respect to analog shaping are that a large number 
of filters are available without restriction to the possible 
settings of the analog shaping module and that raw data 
remain available for further reprocessing. 

2.1.2 Energy Resolution 

The energy resolution of a germanium detector depends 
on the electronic noise, on the charge production in the 
crystal and on the charge collection properties of the 
diode and the shaping filter. A hypothetical 7 line at 
energy E will have a AE (FWHM) expressed by: 

AE = 2.355 \J ^ ENC 2 +r]F-E + c 2 E 2 (4) 

where: 

— 77 is the average energy necessary to generate an 
electron-hole pair ( rj = 2.96 eV in Ge) and F is 
the Fano factor (~ 0.1 for Ge [14]). This term con¬ 
tributes with about 1.8 keV at 2039 keV thus im¬ 
posing a lower limit to the achievable AE; 

— c is a parameter related to the quality of the charge 
collection and integration. An incomplete charge col¬ 
lection can be induced by charge recombination due 
to a too high impurity concentration or due to a too 
low bias voltage while a deficient integration of the 
collected charge can arise if a filter with a too short 
integration time is employed. The same effect is ob¬ 
tained in all cases resulting in low-energy tails of the 
spectral peaks. The parameter c expresses therefore 
the amplitude of such tails. For the detectors used in 
Gerda Phase I, the third term of Eq. 4 is usually 


Table 1 Definition of data sets. The run ranges and active 
detectors are listed. 


set 

duration 

detector configuration 

A 

09.11.11-22.05.12 

ANGs+RGs+GTFs 

B 

02.06.12-15.06.12 

ANGs+RGs+GTF112 

C 

15.06.12-02.07.12 

ANGs+RGs+GTF112 

D 

08.07.12-21.05.13 

ANGs+RGs+GTFl 12+BEGes 


one order of magnitude lower than the electronic 
and charge production terms for events with energy 
up to 3 MeV. 

If the charge collection inefficiency is not dominant, 
the optimization of the energy resolution depends al¬ 
most exclusively on ENC, i.e. on the shaping filter. 
Given that ENC is independent of the energy, any 7 
line with sufficiently high statistics can be exploited for 
the optimization of the shaping filter. 


2.2 Data collection and processing in Gerda 

Calibration data from the period Nov. 2011 - May 2013 
were used to optimize the shaping filter. The detectors 
considered are ANG2-5 from the HdM experiment, 
RG1-2 from IGEX and four of the five BEGes (with 
names starting with “GD”). These are the same de¬ 
tectors used for the 0vj3j3 decay analysis [2]. Since the 
electronic disturbances could change as a function of 
the detector configuration in Gerda, the calibration 
data were divided in four data sets as listed in Table 1. 
In total 72 (45) calibration measurements are available 
for the coaxial (BEGe) detectors. 

2.2.1 Calibration of the energy spectrum 

The calibrations were performed by inserting up to three 
228 Tlr sources in proximity of the detectors [15,16]. The 
total activity of the sources was about 40 kBq at the be¬ 
ginning of Phase I. The duration of the measurements 
was between one and two hours. The energy threshold 
for the calibrations is ~400 keV to reduce disk usage. At 
least ten peaks with energies between 0.5 and 2.6 MeV 
are visible in the recorded spectra (Fig. 4). While all 
peaks are exploited for the calibration of the energy 
scale, only the full energy peaks (FEP) are used in the 
fit of the FWHM as function of energy. This is neces¬ 
sary because the single escape peak (SEP), the double 
escape peak (DEP) and the 511.0 keV line are Doppler 
broadened. 

Given the large number of calibration spectra to be 
analyzed, a fully automatized routine was developed 










5 



Fig. 4 A 22S Th calibration spectrum recorded by ANG5. The threshold is set to ~400 keV. 


and used throughout Phase I. The main steps of the 
procedure are: 

— rejection of events which might decrease the pre¬ 
cision of the calibration; e.g., coincidences between 
detectors, wave forms with superimposed events (pile- 
up events); 

— search and identification of the peaks; 

— fit of the peaks and automatic adjustment of the 
fitting function according to the number of events 
in the peak and the peak shape; 

— extraction of the calibration curve; 

— fit of the FWHM as a function of energy. 

2.2.2 Signal processing 

The signal processing of Gerda Phase I data was per¬ 
formed through an offline analysis of the digitized wave 
forms with the software tool Gelatio [17]. The stan¬ 
dard energy reconstruction algorithm is a digital pseudo- 
Gaussian filter consisting of: 

— a delayed differentiation of the sampled trace 

[t\ -t Xi [t] = x 0 [t] -x 0 [t- 5] (5) 

where Xg [t] is the signal height at time t and S was 
chosen to be 5 ps; 

— the iteration of 25 moving average (MA) operations: 

1 

Xi[t\ -A x i+1 [t] = - ^ Xi[t'] i = 1,..., 25 
° t'—t—5 

( 6 ) 

The energy is given by the height of the output signal 
whose shape is close to a Gaussian. 


This pseudo-Gaussian shaping is a high-pass filter 
followed by n low-pass filters. The resolution obtained 
with the pseudo-Gaussian shaping is very close to op¬ 
timal if the detectors are operated in conditions where 
the 1// noise is negligible [12]. This is not the case 
for Gerda Phase I where the preamplifiers had to be 
placed at a distance of 30-60 cm from the crystals 
due to the low background requirements. The diodes 
and the pre-amplification chain were connected by cop¬ 
per stripes. Hence, a significant low-frequency noise is 
present for some of the Gerda Phase I detectors. 

As described in Sec. 2.1, the ENC depends on the 
properties of the detector, of the preamplifier and of 
the connection between them. In Gerda the diodes 
have different geometries and impurity concentrations 
resulting in different capacitances Cd and different //,. 
In addition, the non-standard connections between the 
detectors and the preamplifiers result in different input 
capacitances (Ci). It is therefore preferable to adapt the 
form and the parameters of the shaping filter to each 
detector separately. 


3 ZAC: a novel filter for enhanced energy 
resolution 

Several methods have been developed to obtain the op¬ 
timum digital shaping for a given experimental setup [12, 
13,18,19]. For series and parallel noise and with in¬ 
finitely long wave forms it can be proven [18] that the 
optimum shaping filter for energy estimation of a 6 - 
like signal is an infinite cusp with the sides of the form 
exp (t/r s ) where r s is the reciprocal of the corner fre¬ 
quency; i.e., the frequency at which the contribution of 
the series and parallel noise of the referred input be¬ 
come equal. When dealing with wave forms of finite 
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length, a modified cusp is obtained in which the two 
sides have the form of a sink- curve. If low-frequency 
noise and disturbances are also present, the energy res¬ 
olution is optimized using filters with total area equal to 
zero [20] . In addition, the low-frequency baseline fluctu¬ 
ations (e.g. due to microphonics) are well subtracted by 
filters with parabolic shape [21], The best energy reso¬ 
lution for Gerda is achieved if a finite-length cusp-like 
filter with zero total area is employed. This can be ob¬ 
tained by subtracting two parabolas from the sides of 
the cusp filter keeping the area under the parabolas 
equal to that underlying the cusp. 

In reality the detector output current is not a pure 
5-function, but has a width of approximately 1 ps. If a 
cusp filter is used, this leads to the effect of a ballistic 
deficit [22,23] and consequently to the presence of low- 
energy tails in the spectral peaks. This can be remedied 
by inserting a flat-top in the central part of the cusp 
with a width equal to almost the maximum length of 
the charge collection in the diode. The resulting filter is 
a Zero-Area finite-length Cusp filter with central flat- 
top that will be referred as ZAC from here on. 

The ZAC filter was implemented as: 


ZACit) = 
sinhf — ) + A ■ 


(*-f ) 2 -(tr 


0 < t < L 


sinhf j- 

sinh ( 2L + FT ~ t i + A ■ 


L < t < L + FT 

(iL + FT-t ) 2 -(f) 

L + FT <t <2L + FT 


( 7 ) 


where t s is the equivalent of the shaping time for an 
analog shaping filter, 2 L is the length of the cusp filter 
and FT is that of the flat-top and where the constant 
A is chosen such that the total integral is zero. The nu¬ 
merical expression of the ZAC filter is obtained through 
the substitution t —¥ At ■ i where At is the sampling 
time and i the sample index; the maximum number of 
samples in the ZAC filter is hzac- A graphical repre¬ 
sentation of the ZAC filter construction is provided in 
Fig. 5. 

Before proceeding with the shaping the original cur¬ 
rent pulse has to be reconstructed from the preamplifier 
output wave form (Fig. 3). This is performed via a de- 
convolution of the preamplifier response function, an 
exponential curve with decay time r = RfCf. Specif¬ 
ically, it is implemented as the convolution with the 


filter consisting of 2 elements, f T 


1 , — exp 




Fig. 5 Amplitude versus time for the ZAC filter (red full 
line). It is composed of the finite-length cusp (blue dashed) 
from which two parabolas are subtracted on the cusp sides 
(green dash-dotted). 


No correction for the finite band width of the electron¬ 
ics was implemented. Since the convolution operation 
is commutative, the convolution between the ZAC filter 
and the inverse preamplifier response function f T can 
be performed once for all: 

FF[i\ = ZAC[i] ■ (-e - ^) + ZAC[i + 1] • 1 

* = 1 , •••) it-zac — 1 ( 8 ) 

The final filter (FF) obtained is shown in red in Fig. 6. 
A convolution of each individual signal trace x with FF 
is then performed: 


i+nzAC —2 

y[i\= ^2 x[k]-FF[i + n Z AC -l~k] 

k—i 

i = 1, ...,n x - n Z AC + 2 (9) 

n x is the number of samples in the trace. Typically n x 
is set to 16384 and uzac ranges from 16060 to 16120. 
The output y for the trace of Fig. 3 is shown as blue 
full line in Fig. 6. The energy E is then estimated as 
the maximum of this convoluted signal y. 


4 Optimization of the ZAC filter on calibration 
data 

The optimization of the ZAC filter using the Phase I 
calibration data was performed separately for each de¬ 
tector. The first and the last calibration run of each 
period were selected (Table 1). Given their longer dura¬ 
tion one more run taken in the middle of the period was 
used for data sets A and D as well. It is expected that 
no change was present in the electronic noise within the 
same data set. In this case the filter parameters giving 
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Fig. 6 The ZAC filter after the convolution with the in¬ 
verse preamplifier response function (red dashed) and the 
wave form of Fig. 3 after the convolution with it (full blue). 


the best energy resolution should be the constant for 
each data set. 

The filter optimization was performed on the FEP of 
208 T1, i.e. the 2614.5 keV line. Quality cuts were applied 
prior to the energy reconstruction that was performed 
only on the surviving events. The energy spectrum was 
reconstructed with different values of the four filter pa¬ 
rameters L , FT, t s and r. In particular: 

— the total filter length 2 L + FT was varied for only 
one calibration run between 120 and 163 ps. As ex¬ 
pected [18] the best energy resolution was obtained 
for the longest possible filter. Given the variability 
of the trigger time within a 2 ps range the maximum 
of the shaped filter can be at one of the extremes 
of the wave form when the maximum filter length 
of 163 ps is used leading to a wrong energy estima¬ 
tion. This effect completely disappears if the filter 
is shortened by 2 ps. Hence, the optimization was 
performed with ~161 ps long filters; 

— the optimal length of FT is related to the charge 
collection time in the detector. For coaxial detectors 
this is typically between 0.6 and 1 ps depending on 
the electric field configuration in the detector and on 
the location of the energy deposition. For BEGes it 
is slightly longer due to the slower charge drift. The 
value of FT was therefore varied between 0.5 and 
1.5 ps in 120 ns steps; 

— the optimal filter shaping time r s depends on the 
electronic noise spectrum as described in Sec. 2.1. 
Typically, r s is of order of 10 ps. The optimization 
was therefore performed with values of r s between 
3 and 30 ps in steps of 1 ps. Since the optimal r s 
was not infinite, the noise present in Phase I data 
had a non negligible parallel component; 

— the value of r can in principle be calculated knowing 
the feedback resistance and capacitance. In reality 
r is modified by the presence of parasitic capaci¬ 


tance in the front-end electronics. Moreover, given 
the presence of long cables a signal deformation can 
arise. Therefore, r is normally estimated by fitting 
the pulse decay tail. This was not possible due to 
the presence of more than one exponential. There¬ 
fore r was varied between 100 and 300 ps with 5 ps 
step size. 


The peak at 2614.5 keV was fitted with the func¬ 
tion [24] for each combination of the filter parameters: 


f{E) = A exp - 


2 a 2 


o C r 

B + — erfc 


(^A 

\V2-a) 


D 

T-exp 


( E- ii 

V 5 


erfc 


E — p 
a/2-ct 



( 10 ) 


corresponding to a Gaussian peak with a low-energy tail 
(last term) sitting on flat background and on a step-like 
function (third term) which describes the continuum on 
the left side of the peak. The FWHM was obtained from 
the fitting function after the subtraction of the flat and 
step-like background components. The energy resolu¬ 
tions resulting from different parameters of the ZAC 
filter were compared and the parameters leading to a 
minimal FWHM were chosen for the full reprocessing 
of the data. For the detectors of the 0pd/3 analysis the 
optimal parameters of the ZAC filter for period D are 
reported in Table 2 as an example. 


5 Results 

The parameter optimization for the ZAC filter provided 
results in agreement with expectations: for each detec¬ 
tor the optimal filter parameters are stable within the 
same data set, but they can vary for those detectors 
that changed configuration in time. This confirms the 
dependence of the microphonic disturbances on the ca¬ 
ble routing. Hence, all Phase I calibration and physics 


Table 2 Optimized parameters of the ZAC filter for pe¬ 
riod D. While the filter length 2 L is equal for all the detec¬ 
tors FT varies between 0.6 and 1.2 ps according to the charge 
collection properties of each diode. 


detector 

2L [ps] 

FT [ns] 

t s [ps] 

T [ps] 

ANG2 

160 

600 

9 

190 

ANG3 

160 

840 

16 

220 

ANG4 

160 

720 

13 

250 

ANG5 

160 

960 

17 

170 

RG1 

160 

720 

12 

210 

RG2 

160 

680 

8 

240 

GD32B 

160 

1080 

13 

220 

GD32C 

160 

960 

16 

170 

GD32D 

160 

840 

15.5 

170 

GD35B 

160 

1200 

17 

135 
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Fig. 7 208 T1 FEP data for ANG5 at 2614.5 keV. The curves 
and parameter values corresponding to the best fit for the 
ZAC and the pseudo-Gaussian shaping are shown. 


Fig. 8 208 T1 FEP data for GD35B at 2614.5 keV. The 

curves and parameter values corresponding to the best fit 
for the ZAC and the pseudo-Gaussian shaping are shown. 


data were reprocessed with the optimal parameters of 
the ZAC filter. 

A first remarkable result is the improvement of the 
energy resolution between 5 and 23 % for the 208 T1 FEP 
at 2614.5 keV of all the Phase I data. As an example 
Figs. 7 and 8 show the summed spectrum of all Phase I 
calibrations around the 2614.5 keV line for ANG5 and 
GD35B, respectively. In both cases, the amplitude of 
the Gaussian component is larger for the spectrum ob¬ 
tained with the optimized ZAC filter and its width is 
correspondingly reduced. The parameters B and C de¬ 
scribing the continuum below the peak are compatible 
for the two shaping filters. 

While for the coaxial ANG5 a low-energy tail has to 
be accounted for in the fit (Fig. 7) the amplitude of the 
tail in the BEGe GD35B is negligible. The tail it there- 

Table 3 Average FWHM over the complete Phase I period. 
The improvement is computed as the difference between the 
FWHM for the pseudo-Gaussian and that for the ZAC filter. 
Only the statistical uncertainty due to the peak fit is quoted. 


detector 

FWHM at 2614.5 keV[keV] 

improvement 

Gaussian 

ZAC 

[keV] 

ANG2 

4.712(3) 

4.314(3) 

0.398(4) 

ANG3 

4.658(3) 

4.390(3) 

0.268(4) 

ANG4 

4.458(3) 

4.151(3) 

0.307(4) 

ANG5 

4.323(3) 

4.022(3) 

0.301(4) 

RG1 

4.595(4) 

4.365(4) 

0.230(6) 

RG2 

5.036(5) 

4.707(4) 

0.329(6) 

GD32B 

2.816(4) 

2.699(3) 

0.117(5) 

GD32C 

2.833(3) 

2.702(3) 

0.131(4) 

GD32D 

2.959(4) 

2.807(3) 

0.152(5) 

GD35B 

3.700(5) 

2.836(3) 

0.864(6) 


fore automatically removed from the fit (Fig. 8). This 
is attributed to the smaller dimensions of the BEGe de¬ 
tector and its reduced charge collection inefficiency. In 
case of ANG5 the tail amplitude D is strongly reduced 
when the ZAC shaping is used thanks to the presence 
of the flat-top that allows for an improved integration 
of the collected charge. 

A deeper understanding of the result is provided 
by studying the evolution of the FWHM as function 
of energy which is fitted according to Eq. 4. An exam¬ 
ple is given in Fig. 9 showing the resolution curve of 



Fig. 9 Resolution curve (eq. 4 with iu 2 = 2.355 2 r/F) for 
ANG5 calculated for all Phase I calibration spectra merged 
together for pseudo-Gaussian (full line) and ZAC (dashed 
line) shaping. 
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Fig. 10 FWHM of the full energy peak of 208 T1 at 2614.5 keV for ANG2. 



Fig. 11 FWHM of the full energy peak of 208 T1 at 2614.5 keV for GD35B. The error bars are partially within the symbols 
size. 


all calibration runs for ANG5. As expected the major 
improvement regards the ENC which reduces FWHM 2 
at all energies by a constant. For both, the pseudo- 
Gaussian and the ZAC filter, the charge production 
term u; 2 = 2.355 2 t?F is compatible with the theoret¬ 
ical value of 1.64 • 1CF 3 keV. Finally, the charge col¬ 
lection term c 2 for the ZAC filter is compatible within 
the uncertainty with the value obtained for the pseudo- 
Gaussian filter. The large uncertainty of this parameter 
is due to the lack of peaks above 3 MeV which makes the 
fit imprecise. This term is the smallest of the three and 
accounts for maximally 15 % of the width at 2614.5 keV. 
A consistent behavior is observed for the other detec¬ 
tors as well. 


One of the original motivations for the application 
of the ZAC filter to the Gerda Phase I data was the 
observation of temporary deterioration of the energy 
resolution in some detectors interpreted as due to time- 
evolving microphonic disturbance not being properly 
treated by the pseudo-Gaussian filter. This is confirmed 
by the comparison of the FWHM over time for both fil¬ 
ters as shown for ANG2 and GD35B in Figs. 10 and 11, 
respectively. In case of ANG2 the FWHM at 2614.5 keV 
obtained with the pseudo-Gaussian shaping fluctuates 
between 4.5 and 4.9 keV. In June 2012 stronger mi¬ 
crophonic disturbance caused a FWHM increase up to 
about 5.1 keV. When using the ZAC filter the effect is 
significantly reduced and the FWHM obtained for the 
affected calibrations is brought back to a value consis- 
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Fig. 12 42 K peak for the coaxial detectors in the energy 

spectrum for all Phase I physics runs. The curves and pa¬ 
rameter values relative to the best fit for the ZAC and the 
pseudo-Gaussian shaping are reported. 


energy [keV] 


Fig. 13 42 K peak for the BEGe detectors in the energy 

spectrum for all Phase I physics runs. The curves and pa¬ 
rameter values relative to the best fit for the ZAC and the 
pseudo-Gaussian shaping are reported. 


tent with the average. Stronger fluctuations were present 
for GD35B. A very poor energy resolution was observed 
during the first month of operation together with a con¬ 
tinuous worsening of the spectroscopic performances in 
the last four months of Phase I. Also in this case the 
ZAC filter energy estimate is unaffected by the low- 
frequency baseline fluctuations induced by microphon¬ 
ics and allowed to stabilize the FWHM over time to 
about 2.8 keV (at 2614.5 keV). 

The Phase I average FWHM for the 208 T1 line at 
2614.5 keV for each detector obtained with the pseudo- 
Gaussian and the ZAC filter are reported in Table 3. 
The average improvement was calculated as the differ¬ 
ence between the two values. This is about 0.31 keV for 
the coaxial and 0.13 keV for the BEGe detectors apart 
from GD35B for which a much larger improvement is 
obtained as described above. 

The comparison of the effective energy resolution 
achieved with Phase I physics data can be performed 
exclusively on the 42 K peak at 1524.6 keV which is 
the only background line with a sufficient number of 
counts for a spectral fit. The summed energy spectra in 
the 1515-1535 keV range for all Phase I data for the 6 
coaxial and the 4 BEGe detectors used for the Oz//3/3 
decay analysis are shown in Fig. 12 and 13, respec¬ 
tively. The FWHM obtained for the pseudo-Gaussian 
shaping and the coaxial detectors is 4.49 ± 0.11 keV. 
This is 0.30 keV larger than the value expected from 
the calibration data. The reason is given by drifts of 
the electronics between calibrations and microphonics 
mainly present in ANG2 and ANG4. For the ZAC filter 


the drifts between different physics runs are reduced be¬ 
cause the microphonics and the noise are treated better. 
The resulting FWHM of the 42 K peak is 4.09±0.11 keV 
and is only 0.15 keV higher than expected from calibra¬ 
tion data. The net improvement in energy resolution at 
1524.6 keV for the coaxial data is 0.40 keV. In case 
of BEGes the ZAC shaping provides a 2.75 ± 0.21 keV 
FWHM compared to 3.05 ±0.30 keV obtained with the 
pseudo-Gaussian. The comparison in this case is harder 
due to the very limited number of events. The improve¬ 
ment on the FWHM of the 42 K line is in agreement 
with the expectation from the calibration data. 

The improvement in energy resolution given by the 
ZAC filter is also reflected in a more precise estimation 
of the energy scale for the single calibration runs. In 
Gerda a second degree polynomial is used as a cali¬ 
bration curve in order to account for the preamplifier 
non-linearity. Figs. 14 and 15 show the residuals of the 
228 Th peak positions from the corresponding calibra¬ 
tion curve averaged over all Phase I calibration runs. 
Both for the Gaussian and the ZAC shaping, the aver¬ 
age residuals are of order of 10” 2 keV. Hence, they are 
much smaller than the peak widths. 

A more informative estimation of the energy cali¬ 
bration precision is obtained by calculating the uncer¬ 
tainty Se of the calibration curve at a given energy, e.g. 
at 1524.6 keV. For each calibration run the quantity 
Se(E = 1524.6 keV) is calculated by error propagation 
on the calibration curve parameters. Using Monte Carlo 
(MC) simulations 10 5 events were randomly generated 
according to a Gaussian distribution with zero mean 
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Fig. 14 Average residuals of the 228 Th peak positions relative to literature values for ANG5. The error bars on the data 
points correspond to the RMS of the residuals for a given peak. 
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Fig. 15 Average residuals of the 228 Th peak positions relative to literature values for GD35B. The error bars on the data 
points correspond to the RMS of the residuals for a given peak. 


and 5e{E = 1524.6 keV). The distributions from all 
Phase I calibration runs are then summed up and the 
systematic uncertainty of the energy scale at 1524.6 keV 
is given by the half-width of the 68 % central interval. 
This results to be between 0.03 and 0.07 keV and is 
up to 16% smaller for ZAC shaping with respect to the 
pseudo-Gaussian filter. 

A cross check of the reprocessed data is given by the 
event-by-event comparison of the energy obtained with 
the ZAC and the pseudo-Gaussian filter. This is per¬ 
formed by calculating the energy difference of the events 
in the 2614.5 keV peak as shown in Fig. 16 for ANG2 
during a typical calibration run. For all the detectors 
this distribution is a Gaussian with a mean value com¬ 
patible with zero and a width <j ~ 0.8 keV. The same 
behavior is observed at all energies for both calibration 
and physics data. 


6 Summary 

The presence of low-frequency noise in the signals of 
Gerda Phase I mostly induced by microphonic distur¬ 
bance leads to a degraded energy resolution for some 
of the deployed detectors. Spectroscopic performance 



Fig. 16 Distribution of the difference between the energy 
estimated with the pseudo-Gaussian and that obtained with 
the ZAC filter for the 208 T1 FEP events at 2614.5 keV. The 
data refer to a standard calibration and are for ANG2. 



























12 


close to optimal is obtained by the use of the ZAC 
shaping filter. This novel Zero Area Cusp filter is ob¬ 
tained by subtracting two parabolas from the sides of 
the cusp filter keeping the area under the parabolas 
equal to that underlying the cusp. A selection of cal¬ 
ibration runs has been exploited for the optimization 
of the ZAC filter. All calibration data sets have then 
been reprocessed using the optimal filter parameters. 

An average improvement of 0.30 keV in FWHM has 
been obtained for both coaxial and BEGe detectors. In 
one case (GD35B) the energy resolution is improved 
by 0.86 keV with the excellent low-frequency rejection 
provided by the ZAC filter. 

The stability of the filter parameters over time for 
the same detector configuration in Gerda along with 
its outstanding low-frequency noise rejection capabili¬ 
ties provides a FWHM improvement of 0.40 (0.30) keV 
at the 42 K line in the Phase I physics data for the coax¬ 
ial (BEGe) detectors. Any improvement in the energy 
resolution will increase the sensitivity of the experiment 
and allow a better understanding of the experimental 
background. 

The Phase I physics data, reprocessed with the ZAC 
shaping, will be combined with the Phase II data in a 
future analysis of the 0v/3(3 decay. The optimization of 
the shaping filter will be performed from the beginning 
of Phase II following a procedure similar to the one 
described in the present work. 
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